COMBINATORIAL ASPECTS OF ORTHOGONAL GROUP 

INTEGRALS 

TEODOR BANICA AND JEAN-MARC SCHLENKER 

Abstract. We study the integrals of type 1(a) = J Q Yi u ij J du, depending on a matrix 
a £ M pxq (N), whose exact computation is an open problem. Our results are as follows: 
(1) an extension of the "elementary expansion" formula from the case a £ M2 X9 (2N) to 
the general case a £ M pxq (N), (2) the construction of the "best algebraic normalization" 
of 1(a), in the case a £ M2 Xq (N), (3) an explicit formula for 1(a), for diagonal matrices 
a £ Af3 X 3(N), (4) a modelling result in the case a £ Mi X 2(N), in relation with the 
Euler-Rodrigues formula. Most proofs use various combinatorial techniques. 



Introduction 

An interesting open question, with several potential applications, is the exact compu- 
tation of the polynomial integrals over the orthogonal group O n . These integrals are best 
introduced in a "rectangular form", as functions of a matrix a G M pxq (N), as follows: 



1(a) = [ fill 1 ' U 11 



Of course, the case p = q = n is the one we are interested in. The parameters p, q are 
there only because they measure, and in a quite efficient way, the "complexity" of the 
problem. As an example, at p — 1 one can use the basic fact that the first slice of O n is 
isomorphic to the sphere S" 1-1 , and the formula for the above integrals is simply 



, v (n-l)!!^!!...^!! 



<?•_ 



where m\\ = (m — l)(m — 3)(m — 5) . . ., stopping at 1 or 2, and where e G {0, 1}. 

The complexity of the problem is measured as well by the value of n G N. For instance 
at n = 2, a classical computation gives the following formula, with e G { — 1, 0, 1}: 

a c\ (a + d)U(b + c)H 



b dl (a + b + c + d+ V 
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Another result, due to Diaconis and Shahshahani [15] . is that the variables n 1 ^ 2 Uij are 
Gaussian and independent at n = oo. This gives the following formula, with e G {0, 1}: 

1(a) = n- k JJa^!! + O^- 1 ) j . 

All these computations, basically done by using geometric methods, will be explained in 
detail in Section 1 below. We refer as well to Section 1 for the precise values of e. Some 
other geometric methods were developed by Olshanskii [19] and Neretin [T7] . 

The integrals 1(a) can be studied as well by using combinatorial methods. It was 
known since Brauer [7] that the intertwining space End(w® fc ) is spanned by certain explicit 
operators, associated to the pairings between 2k points. On the other hand, by basic 
representation theory, the orthogonal projection onto End(ii® fc ) is nothing but the matrix 
formed by the above quantities 1(a), with Ea^ = 2k. So, Brauer's result reduces the 
computation of 1(a) to certain combinatorial questions regarding pairings. The precise 
formula here, whose origins go back to Weingarten's paper [21], is as follows: 

/ v>hh ■ ■ ■ u i2khk du= <*w(i)kO')Wfc»(' r > a). 

Here Dk is the set of pairings of {1, ... , 2k}, the delta symbols are or 1, depending on 
whether the indices match or not, and W^n = G^, where Gf. n (ir, a) = n) 7rS/a '. 

In general, the inversion of Gkn is a difficult linear algebra problem. The first results 
here go back to Collins' paper [8j for the unitary group U n , and then to the paper of 
Collins and Sniady [TT], where several other groups, including O n , are discused. More 
results here were obtained by Collins and Matsumoto in [10], and Zinn- Justin [22]. These 
combinatorial methods have a number of concrete applications, notably to random matrix 
questions [6], [9], [11], [12], [16], [18], [20], and to invariance questions [5]. 

The present paper is a continuation of our previous work [3J. We will use the "elemen- 
tary expansion" approach to the Weingarten function Wk n developed there, in order to 
reach to some new results regarding the integrals 1(a). Our results are as follows: 

(1) We first have some theoretical results extending the "elementary expansion" tech- 
niques in [3], from the case a G M 2X9 (2N) to the general case a G M pxq (N). 

(2) As a first application of our improved elementary expansion methods we will con- 
struct the "best algebraic normalization" of 1(a), in the case a G M 2xg (N). 

(3) In the case where a G M pxp (N) is diagonal the integrals 1(a) were explicitely 
computed in [3] at p — 2. Here we obtain an explicit formula at p — 3. 

(4) Finally, we have a probabilistic modelling result for the pair of variables (itn, ^12), 
in relation with the Euler-Rodrigues formula, which is available at n — 3. 

The precise statements of the above results, as well as a substantial number of techni- 
cal versions and comments regarding what happens in relation with each problem when 
looking at more general p x q shapes will be given in the body of the paper. 
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The paper is organized as follows: Section 1 is preliminary, and in Sections 2-3, 4-5, 
6-7, 8-9 we prove the above results (1-4). Section 10 contains a few concluding remarks. 

Acknowledgements. We would like to thank B. Collins for several useful discussions. 
The work of T.B. was supported by the ANR grants "Galoisint" and "Granma". 

1. Group integrals 

Our basic object of study will be the polynomial integrals over O n . It is convenient to 
introduce these integrals in a "rectangular form" , as follows. 

Definition 1.1. Associated to any matrix a G M pxq (N) is the integral 

v q 



J On ;_1 „ _1 



i=l j=l 

with respect to the uniform measure on the orthogonal group O n . 

Our main tool for the computation of integrals over O n will be a combinatorial formula 
due to Collins and Sniady [TTJ, whose origins go back to Weingarten's paper [2"T] . 

Theorem 1.2. We have the Weingarten formula 

I U ilh . . . U i2khk du= ^ SA^SvUWkniTT, 0"), 

where the objects on the right are as follows: 

(1) D k is the set of pairings of {1, ... , 2k}, also called Brauer diagrams. 

(2) The delta symbols are or 1, depending on whether indices match or not. 

(3) The Weingarten matrix is Wu n = G^ n , where C7fc n (7r,(x) = n^ v ^. 

As an example, the integrals of quantities of type u il j 1 Ui 2 j 2 Ui 3 j 3 u ii j 4 appear as sums of 
coefficients of the Weingarten matrix W^n, which is given by: 

i 2 n n\ -y (n + 1 — 1 — 1 

W 2n = I n n 2 n \ = — — — | -1 n + 1 -1 




n(n-l)(n + 2) 1 _ x _ 1 n + 1 
More precisely, the various integrals at k = 2 can be computed as follows: 

/ u n u n u n u u du = y~] W 2n {^, cr) = , 3 oN , 
Jo n n(n + 2) 

'2 2\ f \ 

n n J = / UuUuU 12 u 12 du = }W 2 n{^,^) = — -, — —77, 
U Uy J Qn ^ n(n + 2) 

2 0\ f n + 1 

u\\u\\u 22 u 22 du = W 2n (r\n, nn) 



2 J J Qn n(n-l)(n + 2) 
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In general, the computation of the Weingarten matrix is a quite subtle combinatorial 
problem and the first results here go back to [21], [TT] . A quite powerful formula was 
recently obtained in [10] and was further processed and clarified in [22J. 

As a first application of the above Weingarten methods let us record a first basic result, 
discussed in [TJ, regarding the vanishing of 1(a), and its n — > oo behavior. We use the 
notation m\\ = [m — l)(m — 3)(m — 5) . . ., with the product stopping at 1 or 2. 

Proposition 1.3. If the sum in some row or column of a is odd, then 1(a) = 0. Also, 



1(a) = n~ k • JJa^l! + 0(n _1 )j , 



where k = Ea^-, and e = 1 if all the entries are even, and e = otherwise. 

We can see from the above result that 1(a) depends in a non-trivial way on the parity 
of the entries a^. We have the following definition. 

Definition 1.4. A matrix a £ M pxg (N) is called "admissible" if the sum on each of its 
rows and columns is an even number. 

In what follows we will restrict attention to the admissible matrices. We should mention 
here that we do not know if the converse of the vanishing statement in Proposition 1.3 
holds or not. This is an interesting question, that we will not investigate in this paper. 

Finally, let us work out a few concrete examples. 

Proposition 1.5. At n = 2 we have the formula 



a c 



i : , = e ■ 



(a + d)\\(b + c)H 
b d)~" ( a + b + c + d+ 1)!!' 

where e = 1 if a,b,c,d are even, e = — 1 is a,b,c,d are odd, and e = otherwise. 

Proof. The idea is that we can restrict the integration to SO2 = S 1 , and the result follows 
from a well-known trigonometric identity. See [2]. □ 

Proposition 1.6. For any a £ N 9 we have 

T . (n-l)\\aA\...a q \\ 
(Edj + n — 1)!! 

where e = 1 if all entries of a are even, and e = otherwise. 



Proof. At n = 2 this follows from Proposition 1.5. In general, the result follows from this 
observation by using spherical coordinates and the Fubini theorem. See [2]. □ 



COMBINATORIAL ASPECTS OF ORTHOGONAL GROUP INTEGRALS 



5 



2. Elementary expansion 



We discuss here a general method from [3] for the computation of the integrals 1(a). 
We will extend this method, from the case a £ M 2xg (2N) discussed in [3] to the general 
case, where a £ M pxq (N) is "admissible" in the sense of Definition 1.4. 

We begin with a number of combinatorial definitions. 



Definition 2.1. Let oi, . . . 

be a pairing of the set {1, . 



, a p £ N be such that the number A 



Safc is even, and let a 
A}. We call "type" of a the matrix r £ M P (N) where rki is 



the number of elements of {(£ 



k-l 



to elements of {(£ 



i-i 



+ i,(Eti«i) + 2,...,(E 



+ i,(£ 



fc-i 



+ 2, • • • , {J2i=i a i) + a k) joined by a 



i-i 



+ at}. 



Observe that this matrix r £ M p (N) is by definition symmetric, even on the diagonal, 
and the sum on its i-th row or column is equal to dj, for any % £ {1, . . . ,p}. 

As an example, let p = 3 and ax = 3, a 2 = 5, 03 = 4. In this case the total sum A = 12 
is indeed even, and the set {1, . . . , A} to be partitioned is: 

{1, . . . , 12} = {1, 2, 3, 4, 5, 6, 7, 8, 9, TO, IT, 12}. 

Here the upper and lower bars are there for reminding us that we have 12 = 3 + 5 + 4. 
Now consider an arbitrary pairing of this set, as for instance the following one: 

o- = { ( 1 , 4) , (2 , 6) , (3, 9) , (5 , 7) , (8, 12) , ( 10, IT) } 

According to the above definition, the type of this pairing is simply the symmetric 3x3 
matrix whose entries count the various types of paired points (underlined-underlined, 
underlined-overlined, and so on). So, the type of this pairing is: 




Now back to the general case, we will need the following observation, which follows 
from definitions. The number of pairings of a given type r £ M p (N) is given by: 



KJa) 



n 

,k<l 



ai 



TkV- 



n 



In what follows we will actually need an abstract formulation plus matrix generalization 
of the above notions. We have the following key definition. 

Definition 2.2. To any column vector a £ M pxl (N) we associate the set [a] consisting 
of the matrices r £ M p (N) which are symmetric, even on the diagonal, and such that the 
sum on the i-th row and on the i-th column is equal to a{. For r £ [a] we let: 



KJa) 



n 

.k<l 



Tkl 



TkV- 



IT 
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More generally, to any matrix a £ M pxq (N) we associate the set [a] = [ai] X . . . x [a q ], 
where a±, . . . , a q are the columns of a, and for r £ [a] we let K r (a) = Ylj=i K r i ( a j) ■ 

We recall from [3] that a matrix a £ M pxq (N) is called elementary when the sum on 
each of its columns is equal to 2. As pointed out there, all the Weingarten function values 
Wkni^-, o~) are quantities of type 1(a), with a being an elementary matrix. Moreover, it was 
shown in [3] that a suitable modification of this idea leads to an "elementary expansion" 
formula for the two-row matrices. We refer the reader to [3J for full details here, the point 
being that a quick reading of [3J is probably useful, but definitely not necessary, for the 
understanding of the most general "elementary expansion" formula, to be explained now. 

For a £ 1R and R £ N we let ot R = a ... a (R times), regarded R vector (that 

we call in the statement below "block"). Also, for r = (r J ) £ [a] we set R^i = E^r^. We 
have then the following result extending the similar "two- row" formula from [3J. 

Proposition 2.3. For any a £ M pxq (N) we have the "elementary expansion" formula 

I(a) = ^2K r (a)I(E r (a)), 

re[a] 

where the elementary matrix E r (a) is constructed as follows: 

(1) There are p rows, and the columns are block-indexed by (k,l) with k <l. 

(2) On the (k,l)-th column (k < I) we have an l Rkl block at rows k and I. 

(3) On the (i,i)-th column we have a 2 R "/ 2 block at row i. 

(4) The other entries are all equal to 0. 

Proof. As a first remark, the above statement might look highly abstract, and its good 
understanding might seem to require a lot of preliminary examples. Indeed, so it is. In 
what follows we will simply present the abstract proof of the above statement, and we 
will leave the examples and particular cases for a bit later. These basically consist of 
the particular case q = 2, discussed in Proposition 3.2 below, and of the particular case 
q = 3, discussed in Proposition 6.1 below. The reader might find it useful to look first at 
Proposition 3.2 and Proposition 6.1 below, as an illustration of the above statement. 
Let us first apply the Weingarten formula. With Ai = X^=i a iji we S e ^ : 

I( a ) = [ <?■■■<? u^.-.u^du 

J On 

= 6 ^ lM ■ ■ ■ P Ap W an ■ ■ ■ Q aiq iapl • • • Q apq W kn (iT, a) 

Tier 

= 5 °( iau ■ ■ ■ 1 aM rpl • • • 1 avq ) 6 ^ lM ■ ■ ■ P Ap ) W kn(K, a). 

(7 7T 

Now let us look at a. In order for the 8 a symbol not to vanish, a must connect between 
themselves the Saa copies of 1, the Ea^ copies of 2, and so on, up to the Ea^ copies of 
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q. So let r G [a] be the type of a in the sense of Definition 2.2. If we denote by [r] the set 
of pairings of type r, the above formula becomes: 

J ( a ) = EEE 6 ^ lAl ■ ■ ■ p Ap ) w u«, <r). 

re [a] <re[r] t 

Let us compute now the elementary integral in the above statement. We label by 01, . . . , <Tk 
the strings of a, and we consider the multi-index j G {1, . . . , k} pk given by j r G o~ r , for 
any r G {l,...,k}, so that 8 a i{j) = 5 ua i for any a'. From the Weingarten formula we get: 

I(E r (a)) = / u lh . . . u ljA u P j Al+ ... +A p _ 1+1 • • • «pja 1+ ...+a p rfM 

= ^^(l A 2 B )Mj)WW7r,a') 

7T(t' 

= ^^(l Al .../-)^ n (vr,a). 

7T 

Together with the above formula and with K r (a) = #[r] this finishes the proof. □ 

As a first consequence of the elementary expansion formula we have the following key 
result, which already appeared in [3] in a more particular formulation. 

Proposition 2.4. We have the "compression formula" 

j fa c\ = Y[cjW J fa £ Cj 



b OJ (£ Cj )H \b 

valid for any vectors a G N q , c G (2N) r and any matrix b G M pxq (N). 

Proof. Consider the elementary expansion formula for the matrix in the statement. By 
splitting apart the indices corresponding to the c columns, we get: 

' (S S) = E * (Z o) ^ (I o c 

re[»] se[g] 

Let us look now at the indices s. These are by definition the r-tuples s = (s 1 , . . . ,s r ) 
with s- 7 G [q] for any i. Now since there is obviously "no combinatorics" when trying to 
pair Cj + + . . . + elements, the only allowed values are s J = diag(cj, 0, . . . , 0), and the 
corresponding coefficients are K Sj Q) = Cj\\. So, with this value of s = (s J ), we get: 

'(S o)-n^E*v(;)/(^.(; „)). 

\ / j rg ra, \ / \ \ / / 
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Let us look now at the matrix on the right. According to the definitions, this matrix is 
obtained from E r (%) by "fattening" it with the Ski = ^j s ii parameters. Now since these 
parameters are simply given by S = diag(Ec.,-, 0, . . . , 0), we conclude that the elementary 
matrix on the right depends only on Ec.,-, so we have a formula of the following type: 

\ / j r6[ a] \ / \ 

We can apply now this formula with Ecj at the place of c, and we get: 

Now by comparing the above two formulae we get the formula in the statement. □ 

The above results suggest the following normalization. 

Definition 2.5. We make the normalization 

(n-gllaj 
1[a) (Eoy + n-l)!!^' 
where I denotes as usual the polynomial integration over O n . 

Let us rewrite the results that we have so far in terms of this new quantity. 

Proposition 2.6. The function J has the following properties: 

(1) Invariance under row and column permutation. 

(2) Invariance under matrix transposition. 

(3) Invariance under horizontal and vertical compression. 

(4) Triviality on the row and column matrices. 

Proof. The assertions (1) and (2) are clear from definitions, and (3), (4) are just reformula- 
tions of Proposition 2.4 and Proposition 1.6, by using the invariance under transposition. 
The "triviality" condition in (4) means of course that we have J G {0, 1}. □ 

We let K' r (a) = K r (a)/([l a i") f° r vectors, and K' r (a) = Yl^r( a j) i n general. 

Theorem 2.7. We have the "compressed elementary expansion" formula 

J(a) = J2K(a)J(E' r (a)), 

re [a] 

where the matrix E' r (a) is constructed as follows: 

(1) There are p rows, and the columns are block-indexed by (k,l) with k < I. 

(2) On the {k,l)-th column (k < I) we have an l Rkl block at rows k and I. 

(3) On the (i,i)-th column we have an Ra entry at row i. 

(4) The other entries are all equal to 0. 
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Proof. This is just a reformulation of Proposition 2.2, by using the compression formula. 
Indeed, consider the elementary expansion formula for the matrix in the above statement: 



1(a) = ^K r {a)I{E r (a)). 



re o 



If we look at the definition of E r (a), we see that for each of these matrices the quantity Ea^- 
is the same as for a, and the quantity Yl a ij * s si m ply equal to 1. Thus when translating 
the above formula in terms of J quantities, the only thing that happens is that the right 
term gets divided by Y[ a^W. This normalization factor can be included into the quantities 
K r (a), and gives rise to modified quantities K' r (a) which appear in the above statement: 

J(a) = J2 K 'M)J(E r (a)). 

rG [a] 

Let us look now at the matrices E r (a) appearing on the right. By using the compression 
formula for each of them, on all rows, we have J(E r (a)) = J(E' r (a)) ) and this finishes the 
proof. □ 

We discuss now some statements which reformulate and generalize the results that are 
left, namely Proposition 1.5 and Proposition 1.6. We use the following notation. 

Definition 2.8. To any aj, b{ G N and n G N we associate the following quantity: 



ai...a k 
6i . . . b s 



(a 1 + n-2)\\...(a k + n-2)\\ 
(6i + n-2)!!...(6 s + n-2)!! ' 



This quantity will be called "balanced" if k = s and Eaj = E6j. 

We already know that we have J = e on one-row vectors, and this reformulates Propo- 
sition 1.6. The following statement is a key generalization of this fact. 



Proposition 2.9. We have the "cross formula" 

/0 6i 
J ci a C2 
\0 b 2 

where B = Efoj, C = Eq 7 and e = 1 if all entries are even, and e 



0(B + C) 
BC 



otherwise. 

Proof. By doing two compressions, it is enough to check the following formula: 



J 



a c 
b 



0(6 + c) 



be 



We proceed by induction, with step 2. At c = 0, 1 the assertion follows from the 1-row 
formula. In general, by using the trivial identity Eufj = 1, we get: 

a c\ _ T fa + 2 c\ f fa c + 2 N 
b ~ [ b + [b 



! ) + <- 2 >^ o o) 
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In terms of J quantities, and after doing a compression, we get: 

( a + 6 + c + n )j^ jj^ =( c + 7i-l)J^ C + 2 ^)+(a + l)J 
Now since by induction the left and right J quantities are equal, we get: 




a c + 2\_b + c + n- l { r a c\ 
b J ~ c + n-1 \b OJ ' 



Thus the formula to be proved is true at c + 2, which finishes the proof. □ 
Finally, Proposition 1.5 has the following reformulation. 

Proposition 2.10. At n = 2 we have the formula 

2 (a + d)(b + c)~ 
abed 

where e = 1 if a,b,c,d are even, e — — 1 if a,b,c,d are odd, and e = otherwise. 

Proof. This is indeed just a reformulation of Proposition 1.5. □ 

3. The two-row case 

In this section we restrict attention to the two-row case, and we extend some technical 
results from [3J, from the "even" case a £ M 2X g(2N) considered there, to the general case 
a £ M 2X(? (N). These extensions will be of great use in Sections 4-5 below. 

We first replace Definition 2.2 by a simpler definition, adapted to the two-row case. 

Definition 3.1. We let [a,b] = {m,m — 2, m — 4, . . .}, with the set stopping at or 1, 
where m = min(a, b). If a, b do not have the same parity, we set [a, b] = 0. 

Observe that this definition agrees with Definition 2.2 in the 2x1 case. Indeed, in both 
situations an element r £ [a, b] is just the type of a pairing of a + b elements. For more 
explanations here, we refer the reader to the proof of Proposition 3.2 below. 

Proposition 3.2. We have the "two-row elementary expansion" formula 

'G)-e*(0'(S vv*). 

^ ' r\...r q \ / \ / 

where the sum is over rj £ [aj, bj], and A = Soj, B = Sfej, R = Erj. 
Proof. We apply Theorem 2.7 to the matrix in the statement: 

'(0-i*G)'(«G))- 



a + 2 
b 
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According to the above new conventions, in this formula the indices range indeed as 
stated, 7j G [dj,bj]. In order to compute the matrix on the right, let us observe that in 
terms of the new indices rj G the matrix indices in Definition 2.2 are: 



„old 



b 3 ~ T l 



By plugging these values into the definition of E' r , we obtain: 

EL 



1 R A-R 
" 1 1 R B — R 



This gives the formula in the statement. 

We denote by A, B,C, . . . the sum of entries of our various matrices a, b, c, 
Theorem 3.3. The function F defined on two-row matrices by 

F (?) = [oiXTB)] J (f) 

satisfies the "flipping formula" F(% — F(f f) for any vectors a, b,c,d. 



□ 



Proof. We apply the elementary expansion formula to the matrices in the statement. Since 
the passage from J to F depends only on the sums on the rows, which are the same for 
the matrix a and for its versions E' r (a), we can replace J by F, and we get: 

/ X R+S A + C _ R _ S 



a c 

b d 

a d 

b c 



EIK 

nsj ij 

EIK 



K' 



K' 



dj 



F ( 



\lR+S 







(l R+s A + D-R 
\l R+s 





B+D-R-S 



S 
B + C-R 



TiSj I] 

Our claim is that the two formulae are in fact identical. Indeed, the first remark is that 
the various indices vary in the same sets. Also, since K' r (£) is symmetric in a, b, the 
coefficients are the same. So, it remains to prove that the F quantities on the right are 
equal. And here, by changing the notations, the statement is that for any s G N and any 
a,b,c,d G N satisfying a + b = c + d, we have: 



I s 2a 
I s 2b 



I s 2c 
I s 2d 



In order to prove this result, we use the cross formula, which gives: 



a c 
b 



(q + c)Q + c) 
c(a + b + c) 
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We prove the result by induction over s. At s = the dependence on a + b follows from 
the cross formula. So, assume that the result is true at s G N. We have: 

V 7 r€[a,6] \ / \ / 

a 0\ v-^ f a\ ^ (\ r a — r 
= 1^ K r 



be ^ r \b \r b + c 

r£[a,b] v 7 

We know from the cross formula that the terms on the left, and hence the sums on the 
right, are equal. With b = s, this equality becomes: 

e k (:) f (;: ■ + 1 - " b °_ r ) = e * (:) p a „ r » + ° - r ) ^ 

r€[o,s] \ / \ / re[o,s] \ / \ / 

Now if we choose a > s, the sums range over the set [a, s] = {s, s — 2, s — 4 . . .}, and by 
the induction assumption, all the terms starting from the second one coincide. Thus, the 
above equality tells us that the first terms (r = s) of the two sums are equal: 

I s a + c-s\ _ (I s a-s 
I s ) ~ \l s c 

Since this equality holds for any a > s and for any c, this shows that the elementary 
flipping formula holds at s, and this finishes the proof. □ 

Proposition 3.4. The quantities in the two-row elementary expansion formula are given 
by 



b) -ll r .! (a ._ r . + i)!!( 6 ._ r . + 

0(A + B-R) 



AB 

where the indices run as usual, rj e [aj,bj], and A = Eaj,.B = Efej,i? = £?y 

Proof. We use the notations from Proposition 3.2 and its proof. According to the identi- 
fications made in the beginning of this section, the coefficients are: 



= n 
= n 



a i !!6^-!(a i -r J -)!r i !(6 J --r i )! 
(a, + !)!!(&, + 1)!! 

r .! (a ._ r . + l)!! (6 ._ r . + l)ir 
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As for the values of J, these can be computed by a method presented in [3]. The idea is 
that, by using several times the flipping principle, these values can be deduced from the 
system of equations coming from the elementary expansion of the quantities J(£) = 1. □ 



Proposition 3.5. We have the "transmutation formula' 

~o(B + cy 



Jn 



a c 
b 



BC 



Jn- 



-C 



valid for any vectors a,b EN q and c G N s . 

Proof. According to Proposition 3.2 and Proposition 3.4, we have: 

'0(A + B + C-R) 



J 71 



Jn 



a c 
b 



+c 



r 



-l) R ' 2 R\\ 
■1) R / 2 R\\ 



(A + C)B 
0(A + B- R)' 



AB 



n+C 



Now observe that the quantities on the right are related by: 



Q(A + B + C - R)~ 




~o(B + cy 




~C(A + B + C - 


-R) 


{A + C)B 


n 


BC 


n 


{A + C){B + C) 






~o(B + cy 




'0(A + B-R)~ 








BC 


n 


AB 


n+C 



Since the other quantities in the above two sums are identical, this gives the result. □ 

Finally, let us record the following useful result concerning some supplementary matrices 
for which / decomposes product of factorials. 

Proposition 3.6. We have the "spark formula" 

"0(a + b + x)(a + b + y) 



J 



x a 
y b 



where e = 1 if a,b,x,y are even, and e 



(a + 6) (a + x)(b + y) 
otherwise. 



Proof. The case e = is clear. In the case e = 1 we use the flipping principle, which tells 
us that the following quantities are equal: 



F 
F 



x a 

y o b 

x 

y a b 



{a + x){b + y) 
0(a + b + x + y) 

x(a + b + y) 
0(a + b + x + y) 



x a 

y o b 

x 

y a b 
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In order to compute the last J quantity we can use the cross formula, and we get: 



J 



x a 0\ xja + b + y) fx 

V b) ~ ( a + x){b + y) \y a b 

0(a + b + x) 
x(a + b) 

After simplifying the fraction, this gives the formula in the statement. 

4. The normalization problem 



x(a + b + y) 
(a + x)(b + y) 

x(a + b + y) 
(a + x)(b + y) 



□ 



We discuss here the normalization problematics for the quantity J(a) for matrices 
a G M 2xq (N) with q = 2, 3. The general case will be discussed in the next section. 
We begin with the case q — 2, where the answer is particularly simple. 



Proposition 4.1. The function tp defined by 

'0 2 {a + d)(b + c) 



J 



abed 



is symmetric in a, b, c, d, and is trivial on crosses, and at n = 2. 

Proof. Since both J and the above normalization constant and invariant under row and 
column permutation, and under transposition, so is the function tp. In addition, we know 
from the flipping principle that we have: 



(a + c)(b + d) 
_0(a + b + c + d) 



J 



(a + d)(b + c) 
_0(a + b + c + d) 



J 



This formula tells us that tp is as well invariant under flliping of columns, so we can 
conclude that tp is symmetric in a, b, c, d. For the cross assertion now, we have: 

"0(6 + c) 



J 



be 



By comparing with the cross formula, we obtain that we have indeed tp 
Finally, the triviality at n = 2 follows from Proposition 2.6. 



e in this case. 

□ 



Our objective now is to find an abstract characterization of the above quantity tp. The 
methods to be developed below will be of great use for the study at q > 3. 

Definition 4.2. A formula of type J(a) = T(a)tp(a) is called "normalization" if: 

(1) r = [Si . . . Sk/Ti . . .T k ], where Si, T; L are sums of entries of a, with Y*Si = STj. 

(2) tp is invariant under column permutation and flip, and is trivial on crosses. 

As a first example, Proposition 4.1 solves the normalization problem at q — 2. 

The above definition probably deserves some more comments. The point is that the 
invariance of a function tp : M2 X q{N) — > M under column permutation and flip, i.e. what 
we can expect from a "careful normalization" of J : M 2X(? (N) — > M, just means invariance 
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under the hyperoctahedral group H q C S^. At q = 2 we have the invariance under the 
whole group S^g, simply because we are in a very special situation, where the transposition 
is available. In the general case q G N, we cannot expect tp to be invariant under S^g- 

We can reformulate the normalization condition in terms of V only. We denote as usual 
by A, B, C, . . . the sum of entries of the various row vectors a,b,c, . . . 

Proposition 4.3. r = [Si . . . S^/Ti . . . is a normalization constant if and only if: 

(1) A(£) = [AB/0(A + B)]T(^) is invariant under column permutation and flip. 

(2) A(l C Q ) = [(a + C)(b + C)/C(a + b + C)}, for any a,be2N and c e (2N) 9 " 1 . 

Proof. We know from Proposition 2.6 that J is invariant under column permutation, so 
the invariance of tp under column permutation is equivalent to that of T, or of A. We also 
know from Theorem 3.3 that [AB/0(A + B)]J(%) is invariant under column flip, so the 
invariance of tp under column flip is equivalent to that of A. Finally, for the checking the 
triviality on crosses we can restrict attention to the case where all the entries are even. 
And here by Proposition 2.9 we have J(% q) = [0(6 + C)/bC], so tp is trivial if and only if: 

0(b + cy 

bC _ ' 

This gives the formula in the statement, and this finishes the proof. □ 

Theorem 4.4. With the above conventions, the normalization problem has a unique so- 
lution for the 2x2 matrices, namely the one given by Proposition 

Proof. Let r(g c d ) = [S± . . . Sk/Ti . . .T k ] be a normalization constant. Here S i: Ti are cer- 
tain sums of entries of c d ). There are 2 4 = 16 such possible sums, and under the action 
of the column permutations and flips, these sums split into 6 classes, as follows: 

(00) = 0, 

(10) = abed, 

(11) = (a + c)(b + d)(a + d)(b + c), 

(20) = (a + b)(c + d), 

(21) = (a + b + c)(a + b + d)(a + c + d)(b + c + d), 

(22) = a + b + c + d. 

Here the symbols on the left, which are all of type (xy) with 2 > x > y > 0, mean that 
we select x entries from the first column, and y entries from the second column. The 
expressions on the right are the products of all sums obtained by using the symbols on 
the left, for the matrix itself, and for its conjugate by column permutation. 
Proposition 4.3 (1) tells us that we must have a formula of the following type: 

A = [(00) A (10) B (11) C (20) D (21) £ (22) F ] . 




(a + C)b 
0(a + b + C) 
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We also know that A must be "balanced" in the sense of Definition 2.8. Thus we have 
the following "balancing equations" , for the total number of terms, and the total sum: 

# : A + AB + AC + 2D + AE + F = 0, 
£ : B + 2C + D + 3E + F = 0. 

Let us look now at the last normalization requirement, namely the triviality on crosses. 
By plugging the value d = into the above expressions, these become: 

(00) = 0, 

(10) = Oabc, 

(11) = ab(a + c)(b + c), 

(20) = c(a + b), 

(21) = (a + b + c)(a + b)(a + c)(b + c), 

(22) = a + b + c. 



According to Proposition 4.3 (2), the triviality on crosses is equivalent to: 

[(00) A (10) B (11) C (20) D (21) S (22) F ] = ^±^±^ . 
By identifying coefficients, we obtain the following supplementary equations: 












A + B = 


= 0, 








a 


B + C = 


= 0, 








c 


B + D-- 


= -1, 






a + 


b 


D + E-- 


= o, 






a + 


c 


C + E = 


= 1, 




a + b + 


c 


E + F = 


= -1. 


unique solution of the system is A = 


C=1,B~- 


= F = 


A = 


[(oo)(ii)- 




~0(a + c)(b + d)(a + 


L(10)(22)_ 




abcd(a + b + c 



■1,D — E — 0, so: 



Thus the normalization constant is: 



r = 



0(a + b + c + d) 
(a + c)(b + d) 



0(a + c)(b + d)(a + d)(b + c) 



2 (a + d)(b + c) 
abed 



abedia + b + c + d) 

But this is exactly the constant in Proposition 4.1, and this finishes the proof. 
We have the following technical extension of Theorem 4.4. 



□ 
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J 



( 



a b c 



) 



abc(a + b + z) (a + c + y) {b + c + x) 
xyz(a + y + z)(b + x + z)(c + x + y) 



( 



a b c 



) 



x y z 



(a + 6) (a + c) (6 + c)(x + y) (x + z)(y + z) 
(a + y) (a + 2) (6 + x) (6 + z) (c + x) (c + y) 



x y z 



is the unique solution of the normalization problem, for the 2x3 matrices. 

Proof. We follow the method in the proof of Theorem 4.4, with less explanations, and 
with some shortcuts. There are 10 symbols of type (xyz) with 2>x>y>z>0, and 
Proposition 4.3 (1) tells us that we must have a formula of the following type: 

A = [(000) A (100) B (110) C7 (111) Z? (200) £ (210) F (211) G (220)^(221) / (222) J ] . 

The "balancing" equations are as follows: 

# : A + 6B + 12C + 8D + 3E + 12F + 12G + 3H + 6I + J = 0, 
E : B + 4C + 4D + E + 6F + 8G + 2H + 51 + J = 0. 

The values of the basic products on a cross (y — z — 0) are: 

(000) = 0, 

(100) = 2 abcx, 

(110) = 0a 2 bcx 2 (a + b)(a + c)(b + c)(b + x)(c + x), 

(111) = ax(a + 6)(a + c)(6 + x)(c + x)(a + 6 + c)(6 + c + x), 
(200) = bc{a + x), 

(210) = bc(a + b)(a + c)(b + c) 2 (a + x) 2 (b + x)(c + x)(a + b + x)(a + c + x), 

(211) = (a + b)(a + c)(a + x)(6 + x)(c + a) (a + 6 + c) 2 

(a + b + x)(a + c + x)(b + c + x) 2 (a + b + c + x), 

(220) = (b + c)(a + b + x)(a + c + x), 

(221) = (a + b + x)(a + c + x)(b + c + x)(a + b + c)(a + b + c + x) 2 , 

(222) = a + b + c + x. 

According to Proposition 4.3 (2), we must have the following formula: 



[(000) A (100) B ...(222) J ] 



(a + b + c)(b + c + x) 
(b + c)(a + b + c + x) 
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We obtain the following supplementary equations: 






A + 2B + C = 0, 


a 


B + 2C + D = 0, 


b 


5 + C + £ + F = 0, 


a + b 


C + £> + F + G = 0, 


b + c 


C + 2F + # = -1, 


a + x 


E + 2F + G = 0, 


a + b + c 


D + 2G + I = 1, 


a + b + x 


F + G + # + / = 0, 


a + b + c + x 


G + 2/ + J= -1. 



The unique solution is A = C = J = 

A = 



-1,B = D = 1,E = F = G = H = I = 0,8q: 
(100)(111) 



(000)(110)(222)_ 
Thus the normalization constant is unique, given by: 

(100)(111) 



"0(^ + 5)" 









(000)(110)(222) 



(iooK(m)Mg) 

(110) 



But this is the constant in the statement, and this finishes the proof. 



□ 



5. Two-row normalization 

In this section we solve the normalization problem for the arbitrary two-row matrices. 
Let us first introduce a candidate for the normalization constant. 

Definition 5.1. To any matrix a G M 2X(J (N) and any r < q we associate the product 



*=nn n z 



Hit 



s seSi s e{i,2} \tes 

with S ranging over all subsets of {1, ... ,q}, having r elements. 

In other words, P r is constructed as follows: (1) we pick r columns of our matrix, (2) 
in each of these columns, we pick one of the two entries, (3) we make the sum of all these 
entries, (4) we make the product of all these sums, over all the possible choices. 



COMBINATORIAL ASPECTS OF ORTHOGONAL GROUP INTEGRALS 



19 



Here are a few examples, for a matrix denoted (£): 
Po = 0, 

! \ <iJ>,- 



Pi 

P2 
P3 



Y\_{di + aj)(bi + bj)(ai + bj)(aj + bA, 



I J (aj + a,j + a k )(ai + a,- + 6 fe )(a, + bj + b k )(bi + bj + 

Observe that all these products are invariant under column permutation and flip. 
Proposition 5.2. For any two-row matrix (£) e M 2xg (N) ; £/ie quantity 

a \ = r_o_ P 9 P 9 - 2 ■ ■ ■ " 

b) [AB' P q -iP q - 3 .. 
where A = Saj,_B = satisfies the balancing assumptions in Definition 2.8. 

Proof. We have to show that for the fraction in the statement, the number of upper terms 
equals the number of lower terms, and the total upper sum equals the total lower sum. 
It is convenient to write our fraction in a closed form, as follows: 




AB 



IP-'- 



r=0 



Since P r is a product of 2 r ( q r ) terms, the signed total number of terms is: 

# = -1 + J2(-iy~ r T ( q ) = -1 + (2 - ly = 0. 

r=0 ^ ' 

Let us look now at the total sums. These are all multiples of the quantity A + B, and for 
r > 1 the averaged total sum inside P r is given by: 



A + B \rj 2q 



SGPr 



r-l 



q-1 
r — 1 



Also, at r = we have S = 0, so the averaged signed total sum is: 

e = -1 + ^(-ir^- 1 (l : j) = -1 + (2 - lr 1 

r=l ^ ' 

Thus T is indeed balanced in the sense of Definition 2.8. 
Definition 5.3. We make the normalization 



0. 



□ 
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where A = Sa^, B = E6j. 

As a first observation, at q = 1 we have the good constant, namely: 

1. 



' 


Pi 




' 


ab 


AB 


Po. 




ab 


~0_ 



At q = 2 now, we recover the normalization constant in Proposition 4.1: 

'0 2 (a 1 + 6 2 )(a 2 + 6i) 



" o P 2 P ~ 




-P*(P 2 /ABY 




AB P x 




Pi 





a x a 2 b x b 2 



Finally, at q = 3 we recover the constant in Proposition 4.5: 



r = 



" o p 3 iV 




[P^Ps/ABY 




Yli a ibi Ui^k^i( a i + a i + b k)( a i + bj + h) 


AB P 2 P _ 




[ Pz 




L Ui^j (a i + a j )(b i + bj ){ai + bj) (a, + h ) 



In what follows we will prove that, in general, tp is indeed a normalization in the sense of 
Definition 4.2, and that it satisfies all the extra properties listed in the introduction. 
In order to prove these results, we use the following key lemma. 

Lemma 5.4. We have the "rational transmutation formula" 

~0(B + C) 



a c 
b 



BC 



r 



a 

" ° 1 b 



valid for any vectors a, b e N q and c eW . 

Proof. We first prove the result at s = 1. Let us introduce the following quantities: 

Pr=Pr(l o), Pr=Pr( a X ^ = U (S + c). 

V J ^ J S€Pr(%) 

Let us look now at the quantity P+. This is the product of sums obtained by selecting 
r columns of the matrix q), then selecting one entry in each of these columns, and 
making the sum. We can see that we have three cases: (1) either the r selected columns 
all belong to the matrix (g), (2) or r — 1 of these columns belong to (£), and the entry 
is selected from the last column, (3) or r — 1 of these columns belong to (£), and the c 
entry is selected from the last column. Thus we have the following formula: 



P+ = P r P r ^P^ v 
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Here we make the convention that both quantities P ?+ i,P_i that can appear are the 
empty products. Now in terms of the T quantity on the left, we get: 



a c 
b 







p + P + 

r q+l r q-l ■ ■ ■ 



(A + c)B P+P+ 2 ... 

(P q P q ^P q ^...)(P^P ( 



( c ) c>(c) \ 
q-2 ■ ■ ■) 



( A + c) B {P q P q -,P q - 2 ...){P^%...) 

n p(c) p(c) 

U r q r q _ 2 . . . 

>(c) 



(A + c)B p( c )p(% 

. v ' ± q—l q — 3 

As for the T quantity on the right, this is given by: 

PqP q -2 



n+c 



AB P 9 _iP 9 _ 3 



(A + c)(B + c) 
Be 



0(B + c)_ 



n+c 

p(c) p(c) 
r 1 r q-2 ■ 

p( c ) p(c) 
r q-l r q-3 

J 



d( c ) p( c ) 
9 - r g-2 



(A + c)P pWpW 



By comparing the above two formulae, this finishes the proof at s = 1. For the general 
case now, s G N, we proceed by recurrence. By using the result at s — 1 and at 1, we get: 

~0(B + C-a) 



a c 
b 



B(C-d) 

o(p + c7- Cl ; 

B(C-ci) 

o(b + c-ci; 

S(C-ci) 



r«+c-ci 

"0(B4 



a Ci 
b 



r 



n+C-ei 



(C-cQCB + CQ 
(P + L7- Cl )C 



n+C 



n+C 



By simplifying we obtain the formula in the statement, and this finishes the proof. □ 

Proposition 5.5. The function (p constructed in Definition 5.3 is a normalization of J, 
in the sense of Definition 4-2. 

Proof. We use Proposition 4.3. The A function is given by: 
A 



(0- 


" AB ' 




' P q P q -2-.. ' 




1 PqPq-2 ■ 




_0(A + B)_ 




[AB P,_!P,_ 3 ..._ 




[A + B P 9 _!P ? _ 3 



Since all the products P r , as well as A + B, are invariant under column permutation and 
flip, so is A. So, it remains to prove that A has the good value on crosses. 
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For a cross (£ g) with a, 5 G 2N and c G (2N) 9 the transmutation formula gives: 





~o(B + cy 




~o(B + cy 




BC 







Thus T has the good value on crosses, and so has A, and this finishes the proof. □ 
We are now in position of stating and proving one of our main results. 

Theorem 5.6. The function ip has the following properties: 

(1) Extension property: (p(%) = (}). 

(2) Invariance under row permutation: </?(£) = 

(3) Invariance under column permutation: tp(%) = ^Cr&j)- 

(4) Invariance under flipping of columns: ip(% = ip(% f). 

(5) Compression property: tp(% q) = tp(% s Cl ). 

(6) Transmutation: (p n (% q) = (p n +c 

(7) Triviality on crosses: tp(% q) G {0, 1} (a,b G N^. 

(8) Triviality atn = 2: (p 2 (t § e {-1,0, 1} (a, b,c,d e N). 

(9) Symmetry at q = 2: tp(% ^) zs symmetric (a, b, c, d G Nj. 

Proof. All the assertions follow from the results that we already have: 
(1) This follows from the rational transmutation formula, at c = 0. 
(2,3,4) These assertions follow from Proposition 5.5. 

(5) The rational transmutation formula shows that T is invariant under compression. 
Since J has as well this property, this gives the result. 

(6) This follows by combining Proposition 3.5 and Lemma 5.4. 

(7) This was already done, in Proposition 5.5. 

(8,9) This was already done, in Proposition 4.1. □ 



6. The three-row case 

In this section we present some advances on the 3x3 case. Let us look first at the 
elementary expansion formula, at p — 3. In this formula the combinatorics comes from 
the pairings of a + b + c elements, and associated to such a pairing are the parameters 
r, s, t, describing the number of "mixed" a — b, a — c, b — c pairings that can appear. 

In terms of the vectors formed by the parameters r, s, t, the formula looks as follows. 

Proposition 6.1. We have the "three-row elementary expansion formula" 

/a\ fa\ fl R I s T A-R-S 

J [b\ =J2Kst\b \ J \ l R s 1 T B-R-T 



rst 



R I s 1 T c-s 



where A, B, C, R, S, T are the sums of entries of the vectors a, b, c, r, s, t. 
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Proof. This simply follows from Theorem 2.7. Observe that the constants appearing in 
the above formula are given by K' rst = Y[K rjSjt:j (aj,bj,Cj)/(aj\\bj\\cj\\), where: 

K rst (a,b,c)=(f\ Q (f) (Y) h (fj r \s\t\(a - r - s)\\(b - r - t)\\(c - s - t)\l 

Indeed, this latter quantity counts the number of pairings of a + b + c elements, having 
r, s, t as the numbers of "mixed" a — b, a — c, b — c pairings. □ 

The above formula shows, and in a very explicit way, that we are in a situation which is 
quite similar to the one in the two-row case, investigated in detail in [3]: the compressed 
elementary matrices are simply the diagonal ones, "decorated" with l's. 

So, in view of the results in [3J, we are led to the following two questions. 

Problem 6.2. What are the values of J on the diagonal 3x3 matrices? And, once these 
values are known, can linear algebra handle the "decoration with 1" operation? 

In what follows we will focus on the first question. In order to deal with the diagonal 
matrices, we use an ad-hoc linear algebra method. Consider the following functions: 

(a 
b 
x V c y 

Here, and in what follows, the matrix entries which are omitted are by definition 0. 
Theorem 6.3. The functions J c (x,y) satisfy the recurrence 

T ( v _ (L + n)J c (x,y)-(x + l)J c (x + 2,y)-(y + l)J c (x,y + 2) 

<J c +2{x,y — : , 

c + n — 2 

where L = a + b + c + x + y, with the following initial data: 

~0(a + b + x)(a + b + y) 



{a + 6) (a + x){b + y) 



Proof. We use the same trick as in the proof of the cross formula. By inserting the trivial 
quantity Ew^ = 1 on the third row, we obtain: 



a \ (a 

b 

x y Cj 
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We perform now the "soft normalization" in Definition 2.5. Most of the constants cancel, 
and we obtain the following formula, where L = a + b + c + x + y: 



(L + n)J 



b 

x y c. 




x y c 2 



The point now is that by doing a compression, we see that the last two J quantities 
appearing on the right are equal. So, the above formula reads: 

(L + n)J c {x, y) = (x + l)J c {x + 2,y) + {y + l)J c (x, y + 2) + {c + n- 2)J c+2 {x, y). 

This gives the recurrence formula in the statement. Finally, the initial data is a value of 
J on a transposed "spark", given by the formula in the statement. □ 

Summarizing, we have a linear algebra method for computing altogether the quantities 
of type J c (x,y), and in particular the "diagonal" quantities of type J c (0, 0). 
As a first application of this method, let us record the following result. 

Proposition 6.4. We have the following formula: 









'1 




1- 









0(a + b) 



(a + 2)(6 + 2) 



1 + 



(a + b + n)(a + n - 2) (6 + n - 2) 
n-2 



Proof. We use Theorem 6.3, with c = x = y = 0. We have L = a + b, and the quantity in 
the statement is therefore given by: 



J 2 (0,0) 



1 



n-2 
1 

n-2 
1 

n-2 



((a + b + n) J o (0, 0) - J (2, 0) - J o (0, 2)) 

0(a + 6 + 2) 



(a + b + n) 
0(a + 6) 



~0(a + b)~ 




ab 





a + b + n) — 



0(a + 6 + 2) 
(a + 2)6 J |_ a(b + 2) 
a+b+n—1 a+b+n—1 



Thus we have a formula of the following type: 

0(a + b) 



J 2 (0,0) 



a + n — 1 



b + n- 1 



(a + 2)(6 + 2) 



n-2 
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More precisely, with S = a + b + n, A = a + n — 2, B = b + n — 2, we have: 

K = 3{A + l){B + l)-(3-l){3 + n-2) 
= S(AB + A + B + 1) - (S- l)(S + n- 2) 
= S{AB + S + n-3) - (5-l)(5 + n-2) 
= SAB + £(£ + n-3)-(£-l)(£ + n-2) 
= SA5 + (n-2). 

This gives the formula in the statement. □ 

The above computation is quite instructive, because it shows that we cannot have a 
flipping principle in the 3x3 case. Indeed, by flipping the "2" entry of the matrix in 
Proposition 6.4 we obtain a cross, and the value of J on this cross is known to be a 
product of factorials. On the other hand, the value of J on the matrix in Proposition 6.4 
is not a product of factorials. Thus, we cannot have a 3 x 3 flipping principle. 



7. Moments of 3 coordinates 



The joint moments of «n,W22 were explicitely computed in [3]. In this section we 
present an explicit formula for the joint moments of Uu, U22, %3- 
We begin our study by solving the recurrence in Theorem 6.3. 



Proposition 7.1. We have the integration formula 



J 



a 

b 

x y c. 



c/2 
k=0 



c/2 
k 



E 

r+s=k 



x + 2r 



X 



2 L 



y + 2s 
y 



j 2 



"0" 




L 




c 


n-1 


_L-c + 2k_ 


n+l 



0(a + b + x + 2r)(a + b + y + 2s) 
(a + b)(a + x + 2r)(b + y + 2s) 

where L = a + b + c + x + y. 

Proof. According to Theorem 6.3, with L = a + b + c + x + y, we have: 

(L + n)J c {x, y)-(x + l)J c {x + 2,y)-(y + l)J c {x, y + 2) 



J c +2(x,y) 

By applying this formula k times, we obtain 

(c-2k + n-l 



c + n-2 



J c+2 (x,y) 



(c + n-1)!! 



E 



fc! 



(L + n + 1)!! 



r+s+t=k 

(x + 2r)H {y + 2s)\\ 



XU 



y\\ 



' r\s\t\ (L-2t + n + l)\\ 
J c -2k+2{x + 2r,y + 2s). 
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With the change c — > c — 2, this formula becomes: 

(c-2fc + n-3)!! 



Jc{x,y) = 



(c + n-3)!! 



E (-D 



fc! 



(L + n- 1)!! 



r+s+t=fc 

(x + 2r)H (y + 2s)!! 



r!s!t! (L-2i + n-l)!! 
^c-2fc(z + 2r,y + 2s). 



X!! 7/!! 

With the choice k = c/2, this formula becomes: 

(n-3)!! ^ , (L + n-1)!! 



E (-D 



(c + n-3)!! 
(x + 2r)H (y + 2s)!! 



X!! y!! 

By using the initial data in Theorem 6.3, we get: 
a 



r\s\t\ (L-2t + n-l)\\ 
J (x + 2r,y + 2s). 



J 



b 

x y c. 



(n-3)!! 
(c + n-3)!! 



r+s+t=c/2 



(c/2)! (L + n-1)!! 
r!s!t! ' (L-2t + n- 1)!! 



(x + 2r)H (y + 2s)H 



X!! 7/!! 

With A; = r + s, this gives the formula in the statement. 
Theorem 7.2. The joint moments 0/^11,^22,^33 « r e (/wen fry 
'a 

fc 



0(a + 6 + x + 2r)(a + 6 + y + 2s) 
(a + 6) (a + x + 2r)(6 + y + 2s) 



J 



k=0 



c/2\ fc! 
2 



fc! 



r+s=fc 



2r 

r 



"0" 




a + b + c 




c 


n-1 


_a + b + 2k_ 


n+l 



□ 



0(a + 6 + 2r)(a + 6 + 2s) 
(a + b)(a + 2r)(b + 2s) 

where J denotes as usual the normalization of I made in Definition 2.4- 

Proof. We use Proposition 7.1 at x = y — 0. The second line of terms is the one in the 
statement, and the coefficient after the first sum is: 





"2r" 




"2s" 


(5 




2 





fc! 



r\s\ 



(2r)!!(2s)!! 



k\ (2r)\ (2s)! _ fc! ^2r 
Hs! ' 2 r r! ' 2 s s! ~~ 2^ 



2s 

s 



Thus we get the formula in the statement. 



□ 



The above formula is of course not very beautiful, but can be quickly implemented on 
a computer. The experiments based on it lead for instance to the following conjecture. 
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Conjecture 7.3. For any a >b > c, the following quantity is an integer: 



"C" 




~b + c 












.0. 


n-1 











n 



We do not know for the moment how to prove this conjecture. Nor do we have an 
advance on the quite similar conjectures formulated in [3], in the 2x2 case. We believe 
that a good framework for the study of these questions is that of the upper triangular 
3x3 matrices, but for the moment we do not have an explicit integration formula here. 



8. Hyperspherical variables 

In this section and in the next one we explore yet another point of view on the compu- 
tation of the integrals 1(a). We will develop here some direct geometric techniques. 

Our starting point is the Euler-Rodrigues formula. This formula parametrizes the 
rotations of M 3 , by points of the sphere 5 3 C I 4 , as follows: 

/x 2 + y 2 - z 2 -t 2 2(yz-xt) 2(xz + yt) \ 

u=l 2(xt + yz) x 2 + z 2 -y 2 -t 2 2(zt - xy) 

\ 2(yt-xz) 2(xy + zt) x 2 + t 2 -y 2 -z 2 ) 

Here u G SO3 is an arbitrary rotation of M 3 , and (x,y,z,t) G S 3 , which is uniquely 
determined up to a sign change, comes from the double cover S 3 = SU2 — > SO3. 

A natural problem, to be explored in what follows, is whether we have similar models 
in higher dimensions. For instance one can ask for a "probabilistic" modelling of the 
standard matrix coordinates of the group SO n+ i with n G N arbitrary, by variables over 
the sphere S" 2 ™ -1 . So, consider the space Mr, with the standard coordinates denoted 
x\, . . . ,Xn- When N is small, or when we will not need all these coordinates, these will 
be simply denoted x,y,z, . . . We denote by S N ^ C M. N the sphere. 

Proposition 8.1. The following variables have the same law: 

(1) Ex 2 - Eyf : S 2n - 1 ->• R, 

(2) x : S n -»■ R. 

Proof. We can write the metric of S 2n ~ l in terms of the metrics of two copies of S" n_1 at 
distance tt/2: 

gzn-i = cos 2 1 ■ g n -i(x) + sin 2 1 ■ g n -i(y) + dt 2 . 
Here t is the "position" between the two copies of S n ~ 1 , which is on one and 7r/2 on the 
other. This shows that the density of t is given by: 

c 2 cos n_1 1 sin n_1 tdt = c sin™" 1 2t dt. 

Here c n is the volume de 5' n_1 , and c is another constant. Now our function is: 

u = cos 2 1 — sin 2 1 = cos 2t. 
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This shows that the law of u is c(l - u 2 ) 11 / 2 ' 1 du, and this gives the result. □ 

The above result is quite interesting, because at n = 2 we can obtain exactly, modulo 
some rotations, the 9 variables appearing in the Euler-Rodrigues formula. This will be 
explained in detail in the next section. For the moment let us present a "complex variable" 
extension of Proposition 8.1, to be used as well in the next section. We denote the 
coordinates of IR 2n by x\, . . . , x n , yi, . . . , y n . We use as well the complex coordinates 
Zi = Xi + iyi, coming from the identification M. 2n = C n . 

Theorem 8.2. The following variables have the same law: 

(1) Y,z 2 : S 2n ~ l ->■ C, 

(2) z : S n ->■ C. 

Proof. Let N = 2n, and be the volume of S N ~ 1 . We will prove that for both variables 
in the statement, if we write them as Z = X + iY, then the law of (X, Y) is given by: 

d^X, Y) = £=^=-(1 - X 2 - F 2 )(«- 3 )/ 2 dXdY. 

2™ C2n 

(1) Let Si, S*2 be the intersections of S 2n ~ l with W 1 , (iR) n . Let gi,g2 be the metrics on 
these spheres, and let 9 be the distance to Si. The metric on S 2 ^ 1 can be written as: 

d9 2 + cos 2 6> • gi + sin 2 6> • g 2 . 

This shows that the law of 9 is given by: 

fi(9) = ^- cos™- 1 9 sin™' 1 9 d9 = % sin™" 1 29 dO. 

We have X = cos 2 9 — sin 2 9 = cos 29, so the law of X is: 

2 2 

p(X) = ° n sin n - 2 29dX = C " (l-X 2 ) (n ^/ 2 dX. 

We now consider X fixed. For z = (zi, . . . ,z n ) G S* 2 ™ -1 , we set m = Re(z)/\\Re(z)\\ and 
w = Im(z) /\\Im(z)\\. Let a be the angle between u and The law of a is: 

c n— 1 . n _2 j 

sin a da. 

Now by definition of Y we have: 

Y = 2 sin cos9{u, v) = sin 26> cos a. 



COMBINATORIAL ASPECTS OF ORTHOGONAL GROUP INTEGRALS 29 
So the law of Y, for X fixed, is: 

tV\ ° n - X SillW ~ 3 a AV 

vxiY) = —7^ dY 

c n sm 29 

c n ~i sin n " 3 a 

■ dY 



c n sin 29 

c n _x (l-r 2 /sin 2 2g)("~ 3 )/ 2 
c n sin 29 

Cn _x (i -x 2 -r 2 )("- 3 )/ 2 



dY 
dY 



C n (1 -X 2 )(™- 3 )/Vl-X 2 
Cn-1 (1-X 2 -F 2 )(- 3 )/ 2 

c n ' (1 -X 2 )(™- 2 )/ 2 

Together with the law of X found above, this gives the result. 

(2) Let (3i be the angle with e\ — (1, 0, . . . , 0), so that x\ = cos j3±. The law of (3\ is: 

/i(/3 1 ) = ^sin™- 1 /3 1 rf/3 1 . 

The law of x\ is therefore: 

= ^sin"- 2 ^ ^ = -^-(1 -rr?)( n - 2 )/ 2 dxi. 

C n +l C n+ i 

Now let xi be fixed, and let f3 2 be the angle with e 2 = (0, 1, 0, . . . , 0). The law of /3 2 is: 

/i, 1 (/3 2 ) = ^ i sin"- 2 /3 2 rf/3 2 . 

We have x 2 = sin /?i cos (3 2 , so the law of x 2 is: 

c n-l . n ~3 o ^2 



d/jL xl (x 2 ) = sin n /3 2 - 



c n sin (3i 

c n -i f l _ ^ 2 \ (n " 3)/2 dx 2 



c n V. sin 2 Pi J sin ^ 

Cn-l (1 - X 2 - X 2 )(™" 3 )/ 2 



2- 



c n (l-x 2 )(^ 2 )/ 2 

Together with the law of x\ found above, this gives the result. □ 

9. The modelling problem 

Let us go back to Proposition 8.1. This gives a model for the S^-hyperspherical law, 
meaning the common law of the standard n + 1 coordinates of S n C M n+1 , as a variable 
over the sphere S 2 " 1 ^ 1 . The point now is that we can obtain a whole family of variables 
over S 2n ~ l which are S^-hyperspherical, simply by using the action of 2n . 
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We recall that the standard coordinates of R 2n are denoted x\, . . . , x n , y± : . . . , y n . So, 
the points of R 2n are the column vectors (*), and the action of an orthogonal matrix 
U G 2n can be written as U(y) = (^), where Ux, Uy are as usual column vectors. 

Proposition 9.1. The following variables have the same law: 

(1) £(U) = EUxf - TAJy} : S 2 ^ 1 -> R, for any U G 2n , 

(2) u i:j : SO n+ i ->■ R, for any i,j G {1, . . . , n + 1}. 

Proof. This is clear from Proposition 8.1, and from the invariance of the uniform measure 
on the sphere S 2n ~ l under the action of 02«, i.e. from d(y) = dU(y). □ 

The 3 diagonal variables in the Euler-Rodrigues formula are obviously all of the form 
£(U), and the same holds in fact for the 6 off-diagonal variables. This key observation, to 
be discussed in detail a bit later, suggests the following general question. 

Problem 9.2. Given a subset K C {1, . . . , n+1} 2 , is it possible to find matrices U 1 ^ G 02n 
with G K , such that the following families have the same joint law: 

(1) Z(U ij ) : S 211 - 1 ->■ R, with G K, 

(2) mj : SO n+1 ->■ R, TOft (i, j) G X. 

As a first remark, depending on the shape of K, we have several problems: "one-row 
problem" , "2x2 problem" , "diagonal problem" , and so on. The (n + 1) x (n+ 1) problem, 
which of course solves all these questions, will be refered to as the "global problem". 

In what follows we will discuss some simple solutions to some of these problems. It is 
convenient to represent a "solution" by the if-shaped array U K = {U tj \(i,j) G K}. 

Proposition 9.3. The lxl problem has the solution U K = (1). More generally, we have 
the solution U K = (U), for any matrix U G 0<m- 

Proof. The first assertion follows from Proposition 8.1, because the variable £(1) is the 
one appearing there. As for the second assertion, this follows from Proposition 9.1. □ 

Let us discuss now the cases n = 1,2, where the modelling problem can be solved by 
using the well-known structure results for the groups 02,03. In order to deal with the 
case n — 1, let us introduce the following matrices: 

These are the 45° counterclockwise and clockwise rotations in the plane. 
Proposition 9.4. The global n = 1 problem has the following solution: 

V K = ( 1 P 



VP"' \ 

More generally, the 16-element subgroup of 2 generated by S 2 and by p produces via 
U — ?■ £(U) the ^-variable correlated set {± cost, ± sin t), repeated 4 times. 
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P 



-1 



_ 1 fx-y 



Proof. Let us denote by x,y the coordinates of M 2 . The actions of the three matrices in 
the statement are as follows: 

'x\ _ fx\ f x \- ^ f x + y 

j))~\y)' p \y) ~ V2\~ x + y y 

Thus the corresponding £ variables are as follows: 

= x 2 -y 2 , 
Z(p) = l/2((x + y) 2 -(x-y) 2 ) 
Zip- 1 ) = l/2((x-y) 2 -(x + y) 2 ) 
So, the solution proposed in the statement is as follows: 

2xy 



2xy, 
-2xy. 



x 2 — y 2 



-2xy x 2 - y 2 



Consider on the other hand the standard parametrization of S0 2 , namely: 



u 



cos t sin t 
- sin t cos t 



We already know fom Proposition 9.3 that all 8 entries of i{U K ) and of u follow a common 
law (the arcsine one). So, we just have to check that the correlation between x 2 — y 2 , 2xy 
is the same as the one between cost, sint, and this is clear. 

Regarding now the last assertion, the first remark is that the group generated by p 
and by the coordinate switch a — (1 J) is the dihedral group D 8 . This group is given by 
D 8 = Z 8 x Z 2 , and its elements are {p k , ap k \k = 0, 1, . . . , 7}. Now by symmetry reasons, 
these 16 elements produce via U — > £(U) the same variables as the 4 elements {1, er, p, p -1 }, 
repeated 4 times. But these latter 4 elements produce the variables {±(x 2 — y 2 ), ±2xy}, 
which are the same as the variables {±cost, ±sint}, and this finishes the proof. □ 

Let us discuss now the case n = 2. We will see that the global modelling problem here 
is more or less equivalent, probabilistically speaking, to the Euler-Rodrigues formula. 
Consider the 45° counterclockwise and clockwise rotations in the -21,-22 planes: 



/ 1 



P = —nz 



1 

V2 



\ 



1 
1 








1 



o\ 



1 



P = —nz 



1 

72 



Let us introduce as well the following matrices: 

/l fl\ 
10 



T = 







1 



(\ 






(\ 

1 


V° 



1 




1 








1 
1 



o\ 


-1 



0\ 

1 
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Observe that these are the two cycles in S3, acting on the last 3 coordinates. 



Proposition 9.5. The global n = 2 problem has a solution of the following type, where 
Pij are various compositions of p, p^ 1 with the permutations in S4: 



U 



K 




Pl3 
P23 



More generally, the subgroup of O4 generated by S4 and by p produces via U — > £(U) a 
number of copies of the variables ±Uij, with € {1,2,3} 2 7 where u e 5*03. 

Proof. According to the definition of r, r _1 , we have the following formulae: 

£(1) = x 2 + y 2 -z 2 -t 2 , 
£(r) = x 2 + z 2 -t 2 -y 2 , 
ar' 1 ) = x 2 + t 2 -y 2 -z 2 . 

The upper matrices in the statement can be taken as follows: 





(0 


1 


1 








(I 





1 








(0 





1 






1 








-1 




1 





1 





1 




1 


1 


-1 








Pl2 = 7! 





1 


-1 





,Pl3 = 


V2 


1 





-1 





,P23 = 










1 


-1 




V 








l ) 






V° 


1 





-v 






V 


1 








The action of these matrices is as follows: 



P12 







(y + z\ 




(A 




fx + 






(A 




(z + t\ 


y 


1 


x — t 




y 


1 


y + 


t 




y 


1 


x-y 


z 




y- z 


Pl3 


z 




x — 


z 


P23 


z 




z - t 


W 




\x + tj 




w 




\y- 






v) 




\x + yj 



Thus the corresponding £ variables are as follows: 

£(p 12 ) = l/2((y + z) 2 + (x-t) 2 -(y-z) 2 -(x + t) 2 ) = 2(yz-xt), 
£(p 13 ) = l/2((x + z) 2 + (y + t) 2 -(x-z) 2 -(y-t) 2 ) = 2(xz + yt), 
£(p 23 ) = l/2((z + t) 2 + (x-y) 2 -(z-t) 2 -(x + y) 2 ) = 2(zt-xy). 

As for the lower matrices in the statement, these can be taken as follows: 

/Oil \ /1 -1 \ /o 

1 

71 



P21 = 



1 
1 

V 







1 









,P31 



V 



1 

1 

V° 1 




1 




1 




,P32 = 



V 



y/2 



1 


V 1 



1 

1 




1 \ 


-1 
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The action of these matrices is as follows: 



P21 



Thus the corresponding £ variables follows: 

£(p 21 ) = l/2((y + z) 2 + (x + t) 2 -(y-z) 2 
£(p 3 i) = l/2((x~z) 2 + (y + t) 2 -(x + z) 2 
£(p 32 ) = l/2((z + t) 2 + (x + y) 2 -(z-t) 2 - 



fx\ 












fx — z\ 




fx\ 


y 


1 


x + t 




y 


1 


y + t 




y 


z 


"71 


y-z 


, P31 


z 


"71 


X + z 


, P32 


z 


W 




\x-t) 




w 




\y-t) 




w 



1 

71 



f z + t\ 

x + 1 

Z - 7 

\x-yj 



[x 



t) 2 ) 

(y-t) 2 ) 
x ~y) 2 ) 



2(xt + yz), 
2(yt - xz), 
2(xy + zt). 



Summing up, the solution proposed in the statement is as follows: 



' x 2 + y 2 - z 2 - 
2(xt + yz) 
2{yt — xz) 



t 2 2{yz - xt) 
x 2 — y 2 + z 2 — t 2 



2(xz + yt) 
2{zt — xy) 
x 2 — y 2 — z 2 + t 2 



2(xy + zt) 

But this is exactly the Euler-Rodrigues formula, and this finishes the proof. □ 

Theorem 9.6. The 1x2 modelling problem has the solution U K — (l p) , where p is 
the 45° rotation in each Z{ plane. 

Proof. This is just a reformulation of Theorem 8.2, by taking real and imaginary parts, and 
by using the various notations and identifications from the beginning of this section. □ 

10. Concluding remarks 

We have seen in this paper that the computation of the integrals 1(a) leads to a number 
of concrete questions regarding the symmetries and normalization of such integrals, as well 
as to a number of questions regarding their probabilitic modelling, over spheres. 

All these questions depend in a crucial way on the shape p x q of our matrix a. 

Most of the open problems that we have at this time concern the case of 3 x 3 matrices. 
So, the main problem that we would like to raise here is that of understanding the integrals 
1(a) for matrices of type a G M 3x3 (N), perhaps taken to be upper triangular. 

In addition, we have the question of finding noncommutative analogues of these results 
and techniques. The subject here is quite interesting as well, because it is closely related 
to the meander determinant problematics, considered in [13], [TJ]. See [I]. 

We intend to come back to these questions in some future work. 
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